Typeset with jpsj3.cls <ver.l.O> 



Full Paper 



o 
< 



3 



Effect of Disorder on Fermi surface in Heavy Electron Systems 

Junya Otsuki*, Hiroaki Kusunose 1 , and Yoshio KURAMOTO 

Department of Physics, Tohoku University, Sendai 980-8578 
1 Department of Physics, Ehime University, Matsuyama 790-8577 

(Received August 3, 2010) 

The Kondo lattice model with substitutional disorder is studied with attention to the size of 
the Fermi surface and the associated Dingle temperature. The model serves for understanding 
heavy-fermion Ce compounds alloyed with La according to substitution Ce^Lai-a;. The Fermi 
surface is identified from the steepest change of the momentum distribution of conduction elec- 
trons, and is derived at low enough temperature by the dynamical mean-field theory (DMFT) 
combined with the coherent potential approximation (CPA). The Fermi surface without mag- 
netic field increases in size with decreasing x from x = 1 (Ce end), and disappears at such 
x that gives the same number of localized spins as that of conduction electrons. From the 
opposite limit of x = (La end), the Fermi surface broadens quickly as x increases, but stays 
at the same position as that of the La end. With increasing magnetic field, a metamagnetic 
transition occurs, and the Fermi surface above the critical field changes continuously across the 
whole range of x. The Dingle temperature takes a maximum around x = 0.5. Implication of 
the results to experimental observation is discussed. 

KEYWORDS: Large Fermi surface, Dingle temperature, metamagnetism, Ce a; Lai_ :E Ru2Si2, coherent po- 
tential approximation (CPA), continuous-time quantum Monte Carlo (CT-QMC) 
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1. Introduction 

The Fermi surface (FS) reflects strongly the forma- 
tion of heavy quasiparticles. Suppose that 4f electrons 
in Ce compounds behave as localized spins decoupled 
with conduction electrons, the volume surrounded by 
the FS does not include 4f electrons in a sense of Lut- 
tinger's theorem. ^ In this case, the topology of the FS is 
the same with that of the corresponding La compounds, 
which is referred to as "small FS". On the other hand, 
in the heavy-fermion state, the quasiparticles involve 4f 
electrons, which therefore contribute to the FS volume 
to yield "large FS". Thus, investigating the FS is a di- 
rect way of observing the formation of the heavy-fermion 
state. 

CcRu2Si2 is a typical material that exhibits heavy- 
fermion behavior. Magnetic field gives rises to a meta- 
magnetic crossover 2 ) and simultaneously, the topology 
of the FS changes from 4f-itincrant one to that simi- 
lar to LaRu2Si2, accompanied by an enhancement of the 
cyclotron effective mass. 3 -* According to Luttinger's the- 
orem, the total volume surrounded by the FS cannot 
be changed continuously in an ideal system. Hence, the 
change of the FS topology should be referred to as "tran- 
sition" , but should be distinguished from the ordinary 
phase transition in thermodynamics. A similar FS tran- 
sition has been observed by applying pressure in some 
compounds such as CeRhlns 4 ) and Cel^. 5 ) 

Substitution of Ce ions with La ions is an effective 
technique to reveal the role of 4f electrons in physical 
properties. Especially for the FS, it is useful since the 
topologies at both limits are different. Strictly speaking, 
the momentum distribution function has no discontinuity 
even at T = in the alloys with disorder. Nevertheless, 
de Haas-van Alphen (dHvA) oscillations have been ob- 



served for wide range of the Ce concentration. 6 ~ 9 ' Hence, 
in this paper, we use the term "Fermi surface" for the 
steep slope in the momentum distribution. 

Theoretical studies have been performed for binary al- 
loys of Ce and La ions, and the change from the Kondo 
singlet in the dilute limit of Ce to the heavy-fermion state 
in the high density limit has been discussed. 10 ~ 2 °) How- 
ever, the location of the FS in the disordered system has 
not been derived despite extensive studies from exper- 
imental side. In this paper, we theoretically investigate 
the effect of disorder on the FS. In considering dHvA 
results, we should pay attention that dHvA experiments 
observe electronic states under magnetic field. Actually, 
under magnetic field, the FS topology may look as if no 
4f electrons contribute, even though the quasiparticle is 
composed mainly of 4f electrons. 21, 22 ) We shall demon- 
strate that the substitution leads to qualitatively differ- 
ent FS depending on strength of the magnetic field. 

We examine the effect of disorder based on the Kondo 
lattice model. Details of the model and the computa- 
tional method arc explained in the next section. We then 
show numerical results for the case of Ce end (x = 1) 
under magnetic field in §3. Results for disordered system 
are given in §4. We summarize in §5 with discussions on 
realistic materials. 

2. Model and Computational Method 

We consider a disordered Kondo lattice model in which 
localized spins are randomly distributed with concentra- 
tion x. We refer these sites to the A-sites. The Hamilto- 
nian then reads 



*E-mail address: otsuki@cmpt.phys.tohoku.ac.jp 



2 J. Phys. Soc. Jpn. 



Full Paper 



J. Otsuki et al. 



E 



J J2 C itr ffC * 



(Si) e 



2HSf 



(1) 



where a = ±, and c; CT = N~ x l 2 ^ fc Ckae lk Ri with TV be- 
ing the number of sites. The chemical potential includes 
the Zeeman splitting as = fj, — aH, and the g-factor 
for conduction electrons is assumed to be the same as 
localized spins. 

A randomly disordered alloy can be dealt with by the 
coherent potential approximation (CPA). 23,24 ) On the 
other hand, strong local correlation can be taken into 
account by the dynamical mean-field theory (DMFT), 25 ' 
which can be combined with the CPA. 10-19 ) We shall give 
a brief review of this framework for application to the 
Kondo lattice model. 

We assume, for the electronic state, the translational 
symmetry of the lattice. This assumption can be real- 
ized by taking an average over the random configuration 
of two ions. Within a local approximation for the self- 
energy, the single- particle Green function for conduction 
electrons is given by 



Gv(fe, z) = [z - e k + \i a - T, a (z)}~ 



(2) 



The self-energy T, a (z) incorporates both the coherent po- 
tential and the many-body effect due to local interac- 
tions. We find an optimal self-energy from the solution 
of a single-impurity problem in an effective medium. We 
define a cavity Green function Go<r( z ) by 



GoAz)- 1 =G CT (z)- 1 + E CT (z), 



(3) 



where G{z) = TV" 1 J^ fc G(fc, z). We omit the argument 
z and index a hereafter for simplicity. The local Green 
function Q a with a — A, B is defined for the local mo- 
ment site A and the non- magnetic site B. It is given by 



Go + GotaGo, 



(4) 



where t a is the single-site t-matrix, and the interaction 
has been assumed to be local. 

In the CPA, the cavity Green function Go is determined 
so that the local Green function, (G), averaged over two 
kinds of sites is identical to that of the whole lattice, G: 



G={S)= xQa + (1 - x)Ge 



(5) 



This CPA condition can be written in a form which does 
not include Go- Using eq. (3), eq. (4) is rewritten in terms 
of G as 



where 



"a 



G a = G + Gt a G, 

(gg - g) 
l-G(E a -E)' 



(6) 



(7) 



with E" 1 = t,- 1 + Go- The CPA condition (5) is then 
written as 



(t) = xt A + (1 - x)i B = 0, 



(8) 



which leads to 



E = 



(t) 



(t) =xt A + {l-x)t B . (9) 



l + £o<t>' 

This equation gives the self-energy in the CPA. For x 



1, the equation is reduced to that in the pure Kondo 
lattice model. 26 - 1 The effective medium Go is determined 
self-consistcntly from eqs. (2), (3) and (9). 

We solve the effective impurity problem using the 
continuous-time quantum Monte Carlo (CT-QMC) 
method adapted to the Kondo model. 27 ' We use a semi- 
circle density of states p(e) = (2/-7T-D)yl — (e/D) 2 with 
the band cutoff D = 1, and use the following numerical 
parameters: n c = 0.9 for the number n c of conduction 
electrons per site, and J = 0.34. The resultant Kondo 
temperature is estimated as Tk ~ 0.1. With these pa- 
rameters, the antifcrromagnetic ordering is suppressed 
by formation of the Kondo singlet. 28 ) 

3. Effect of Magnetic Field at x = 1 

We begin with the pure system (Ce end) under mag- 
netic field. Figure 1(a) shows the magnetization process. 
We define Hi and H 2 as lower and upper fields of the 
magnetization plateau, respectively The appearance of 
the magnetization plateau is explained in a quasiparticlc 
band picture, which also gives an account of the change 
of the FS topology 21,22 - 1 In the hcavy-fermion system, an 
almost dispcrsionless quasiparticlc band is located close 
to the Fermi level. Mixing with a conduction band re- 
sults in a hybridization gap, which is of the order of Tk 
in the Kondo lattice. 29 ' Under the magnetic field, one 
of spin components is filled up to the energy gap. The 
magnetization does not change in this region, namely, 
H 2 — Hi ~ Tk- 

Let us observe the change of the FS topology. In the 
DMFT, fc-dependence of the Green function enters only 
through £fc as shown in eq. (2). Hence, we regard k = e k 
as "momentum" hereafter. The Fermi momentum can 
be traced by the renormalized chemical potential /i a — 
ReE CT (0) , since the FS of each spin component appears at 
the momentum where n = — RcE CT (0) is satisfied. 29 ' 
Figure 1(b) shows ReE CT (0) as a function of H . There 
are three regions, and correspondingly we introduce three 
"momenta" : kl , us an d n sc , which give a FS involving 
(n c + l)/2, n c /2, and n c electrons, respectively. 
Region I (0 < H < Hi): Without magnetic field, [i a — 
ReE CT (0) corresponds to kl and the FS involves localized 
spins. Upon turning on the magnetic field, the FS splits, 
but the average keeps at «l- 

Region II {Hi < H < H 2 ): At Hi ~ 0.006, one band 
passes the band edge k c = 1 indicating that the chemical 
potential is located in the energy gap. Then, only the 
other spin forms the FS, which has the momentum k sc . 
Namely, the FS of a single spin-component involves n c 
electrons in this region, and the magnetization shows the 
plateau. 

Region III (H > H2): By further increasing the magnetic 
field, both spin components again contribute to the FS. 
For H > H2 — 0.126, the spin average of — RcE (T (0) is 
located at Kg , indicating the small FS composed of both 
spin-components . 

We note that even in Region III, the quasiparticlc 
involves the localized spins near H2, since the energy 
gap cannot be formed without coupling with localized 
spins in the Kondo lattice. Namely, the "small FS" 
does not mean the localized character of 4f electrons. 22 ' 
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the renormalization factor z„ defined by 



Fig. 1. (Color online) Magnetic field dependences of (a) magneti- 
zation M, (b) Fermi momentum re = fi^ — ReS CT (0), and (c) spin- 
dependent renormalization factor z a . The Fermi level lies in the 
energy gap of the majority spin band for Hi < H < H2- Param- 
eters are chosen as x = 1, J = 0.34, n c = 0.9, and T = 0.0025. 



The energy gap gradually vanishes with decreasing 4f- 
contribution to the quasiparticlc for H ^> H2 ~ 7k. 

Transition of the FS topology is known as Lifshitz tran- 
sition, which is distinguished from thermodynamic phase 
transition. 30,31 ) The transitions between Regions I and 
II, or II and III in the Kondo lattice are identified as the 
Lifshitz transition in the sense that there is no disconti- 
nuity in physical quantities at finite temperature. How- 
ever, as important characteristics in the Kondo lattice, 
the transition is caused by quasiparticles, which incorpo- 
rate strong local correlation. This difference appears in 



1 



cTmEo^ie^ 



(10) 



inverse of which corresponds to the mass enhancement 
of the quasiparticle. Since now the energy gap is due to 
the collective Kondo singlet, ImS(w) has a peak centered 
at the energy gap as in the Mott insulator. 25 ' This peak 
leads to a steep slope in the real part ReS(w). As a re- 
sult, the renormalization factor takes a small value when 
the quasiparticle level crosses the Fermi level, as shown 
in Fig. 1(c). Consequently, the effective mass for the ma- 
jority spin is enhanced toward the magnetic field where 
the FS transition occurs, while the other spin, which is 
not involved in the FS transition, becomes lighter on 
the plateau region. The effective mass strongly depends 
on the spin component on both sides of the transition, 
while the spin dependence is negligible for H 3> i?2- 
This enhancement has been discussed in the high-field 
region. 32 ~ 34 ' 

4. Effect of Disorder 

Now let us examine the effect of disorder on the FS. We 
discuss the momentum distribution of conduction elec- 
trons 



( C L- C *^) = T^G cr (K,ie„)c I 



ie„0+ 



n crT (n). (11) 



Figure 2 shows n c+ (n) = n c _(«;) = n c {n) for H = 0. 
Concerning < x < 0.9, the temperature is much lower 
than Tk, and the results may be regarded as those in the 
ground state. In the case of x = 0, n c (n) has a discontinu- 
ity of unit magnitude according to the Fermi-Dirac dis- 
tribution function. The Fermi momentum is determined 
only by the number of conduction electrons. At x = 1, 
on the other hand, a small discontinuity appears at the 
momentum where the Fermi volume includes localized 
spins. 29 ' This large FS disappears at x = 0.9 = n c . The 
volume of the large FS increases by removing localized 
spins toward x = n c . This trend is in contrast to a rigid 
band picture that the volume decreases in proportion 
to the number of localized spins. It can be understood 
by considering that x = n c is a condition for insulator. 
Namely, the Fermi momentum varies toward k c = 1 at 
x = n c . 

For magnetic field of Hi < H < H2, n ca (n) depends 
on n in a similar manner as that in Fig. 2, except that the 
FS in the Ce side appears only in the minority spin. On 
the other hand, when H ^> H2, the FS can be connected 
continuously from x = to x = 1. Figure 3 shows n ccr (K) 
for H = 0.2, which is much larger than H2 and hence 
much larger than Tk- The FS's of majority and minority 
spins split as the Ce concentration increases. The size of 
each component is opposite to the corresponding number 
of local spins because of the effective field through the 
exchange interaction. 

We define an energy scale T ia that characterizes the 
scattering rate of the quasiparticles: 29 ' 

-1 



T, 



(12) 



FS( 7 ) 
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Fig. 2. (Color online) Momentum distribution for H = 0. The 
temperature is taken as T = 0.001 for 0.9 < x < 1 and T = 
0.0025 for < x < 0.88, which are low enough to neglect the 
broadening by finite temperature in the scale of the figure. The 
(red) squares on the bottom and in the inset show the momenta 
where — dn c (re)/dK takes a local maximum. 




Fig. 3. (Color online) Momentum distribution n crT (K) for H = 0.2 
and T = 0.0025. Upper and lower panels correspond to majority 
and minority spins, respectively. The (red) squares on the bottom 
show the momentum where — dn C(T (K)/dK takes the maximum. 

where 7 = S,L, and the subscript FS(7) indicates that 
the derivative is evaluated at the corresponding FS. At 
x = 0, the FS gets blurred according to Ts a = T, since 
n c («) is given by the Fermi-Dirac distribution function. 
On the other hand, the FS at x = 1 depends more 
strongly on temperature: T^ a = T/z 2 . 29 ) Except for 
x = and x = 1, T ia is finite even at T = because of 
disorder. In this case, T 7<T can be regarded as the Dingle 
temperature. 




.V 



Fig. 4. (Color online) Dingle temperature T 7 as a function of x 
for several values of H. For finite H, the open and closed symbols 
represent majority and minority spins, respectively. 

Figure 4 shows T ia as a function of x for several val- 
ues of H . At H = 0, Tg monotonously increases as x 
increases. Namely, the La FS is not connected to the Ce 
side, although it leaves a trace of the "small FS" . Instead, 
the "large FS" appears at x > n c , appearance of which 
is indicated by rapid decrease of Tl. On the other hand, 
when H > H2, To, becomes zero at both x = and x = 1 
having a maximum at 0.5 < x < 0.8. Namely, the FS in 
the La side remains a true FS up to the Ce end. 

5. Summary and Discussions 

We summarize the above results in Fig. 5. The Ce side 
is classified into three regions according to the FS topol- 
ogy: (I) large FS, (II) single-component FS, and (III) 
"small" FS. This transition is caused by the quasiparti- 
cle band, whose energy gap H2 — Hi is of the order of 
Ik. The effect of disorder on the FS is qualitatively dif- 
ferent depending on whether H is larger or smaller than 
i?2- When H < H2, the heavy- fermion FS of one spin 
component appears around x = n c upon increasing x. 
In addition, if H < Hi, the other spin component shows 
up for larger x. In the special case of H = 0, both spin 
components appear simultaneously at x = n Cl at which 
the system becomes insulating. In this magnetic-field re- 
gion, H < H2, the Dingle temperature of the La FS 
monotonously increases against x. On the other hand, 
when H 3> J?2, the La FS is connected continuously to 
the Ce side. In this case, the Dingle temperature has a 
maximum at a certain value of a;. As shown in Fig 5, 
the maximum point tends to converge at x = 0.5 in high 
field, where the system can be regarded as a nonmagnetic 
binary alloy. 

In realistic situations, the chemical potential is not in 
the hybridization gap, because of an energy dispersion of 
4f level and/or existence of plural conduction bands on 
the Fermi level. In this case, a situation may occur where 
the magnetization-plateau region disappears in Fig. 5. 
Then, we expect a metamagnetic behavior without the 
plateau around Hi ~ H2 ~ Ik- 

We finally discuss realistic materials such as 
Ce a; Lai_ a ;Ru2Si2. We consider that the FS transition in 
pure CeRu2Si2 is the Lifshitz transition caused by the 
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Fig. 5. (Color online) The FS topology in x-H plane for J = 0.34 
and n c = 0.9. The open square (red line) represents H of mcta- 
magnctic behavior. The closed square (blue line) represents x 
above which the indication of the quasiparticle is seen in n co -(re). 
The dashed line and circle symbol represent x where Ts CT takes 
the maximum. The open (red) and closed (blue) symbols signify 
anomaly due to majority and minority spins, respectively. 



quasiparticle band. 22 - 1 Actually, an analysis of transport 
coefficients indicates disappearance of a FS with one of 
spin components under magnetic field. 35 - 1 Our numerical 
calculation demonstrates that this FS transition gives 
rise to a metamagnetic transition as shown in Fig. f (a). 
In considering substitution effect in the multi-band sys- 
tem, the Dingle temperature is more informative rather 
than the dHvA frequencies, since each frequency may 
vary with keeping the total volume of the FS. Recent ex- 
periments have revealed a contrastivc behavior depend- 
ing on direction of the magnetic field: dHvA signals van- 
ish by the substitution on both La and Ce sides for 
H _L c, while the dHvA amplitude A osc shows a min- 
imum at < x < 1 for H || c. 36 ' Because the magne- 
tization in this system shows a strong Ising anisotropy, 
we can regard the results for H _L c as low-field region, 
and the results for H || c as high- field region. Then, 
the tendency is consistent with our results for the Dingle 
temperature T 7 , considering A osc oc exp(— T 7 ). The de- 
gree of x-dcpendcncc in the Dingle temperature indicates 
how much 4f electrons contribute to each branch. 
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